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Abstract 

We construct the generating functional of two flavor chiral perturbation theory 
including the effects of virtual photons in the one loop approximation. As an ap- 
plication, we calculate the electromagnetic corrections to the elastic tttt scattering 
amplitude, in particular to its S-wave threshold parameters. Numerical estimates 
are given for the reaction 7r°7r° — > 7r°7r°. These electromagnetic effects are found to 
be smaller than the hadronic two-loop corrections for the scattering length ao- The 
effective range bo increases by 36% due to the unitary cusp. 
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1. Elastic pion-pion scattering is the purest reaction to test our understanding of the 
spontaneous and the explicit chiral symmetry breaking in QCD. This reaction can be 
calculated to high precision within the framework of chiral perturbation theory (CHPT) 
and is one tool to pin down the size of the scalar quark condensate B = \(0\qq\0)\/F , 
with Fq the pion decay constant (in the chiral limit), and to give bounds on the up 
and down quark masses (for a review, see CHPT is the effective field theory of the 
standard model at low energies, which admits an expansion in small external momenta and 
quark (pion) masses, with q collectively denoting any one of these small parameters. This 
so-called chiral expansion, which proceeds in steps of powers of q 2 , can be mapped one- 
to-one on an expansion in the number of pion loops and higher dimension local operators 
0. At present, elastic tttt scattering in the threshold region has been investigated at tree 
level (q 2 ) @], to one loop (g 4 ) [[§] and to two loop accuracy (g 6 ) 0. There exits also a 
dispersive analysis to order g 6 in the framework of generalized CHPT 0. Of particular 
interest are the S-wave scattering lengths a^, with / = 0, 2 the total isospin of the two- 
pion system, since they vanish in the chiral limit 0. Consider as an example a(j in CHPT: 
For the central values of the various input parameters, the tree, one- and two-loop results 
are a{j = 0.156, 0.201 and 0.217, in order .p 5 ! While the one-loop corrections are sizeable 
(~ 25%), the two-loop corrections are already considerably smaller (~ 10%). Apart 
from these strong interaction corrections, there are also electromagnetic contributions. In 
particular, the charged to neutral pion mass difference, which is an effect of the order of 
(M,r+ — M n o)/M n + ~4%, is almost entirely of electromagnetic origin [|7||. It is therefore 
of importance to get a handle on the electromagnetic corrections for the tttt scattering 
amplitude to further sharpen the theoretical predictions. Virtual photons have already 
been included in three flavor CHPT ||0 and various effects like e.g. the violation of 
Dashen's theorem or the electromagnetic (em) corrections to the decay 77 — ^ 37r [K| have 
been calculated. Since Tin scattering can be described purely within SU(2) CHPT, we 
construct here the generating functional including virtual photons to one loop for the 
two flavor case. As proposed in ||, we assign the chiral dimension one to the electric 
charge, so the lowest order effective Lagrangian coupling the virtual photons to the pions 
is of 0(e 2 ) ~ 0(q 2 ). The difference to the three flavor case considered in refs. ||fj 
comes largely from the fact that for two flavors there are less independent operators. 
This reduced number of terms is the natural basis to consider purely pionic processes. 
As an example, we then estimate the em corrections to the tttt scattering amplitude, in 
particular to the threshold parameters (scattering lengths, effective ranges). In this letter, 
we concentrate on the reaction 7r°7r° — > tt°tt since space forbids a thorough discussion of 
the treatment of the infrared divergences appearing in charged pion scattering. Previous 
work mostly related to the leading em corrections involving at least one tt + tt~ pair can 
be found in ([□]] [|T^] [|T3| [|nj . In most of these calculations, the pion mass difference has 
not been included systematically. 

2. The effective field theory build of pions, collected in U(x) = u 2 (x), photons (A^) and 
other scalar (s), pseudoscalar (p), vector (v^) and axial-vector (a M ) external sources starts 

# 5 In particular, = 93.2 MeV and the charged pion mass are used. For a discussion on the parameter 
sensitivity, see H. 



2 



at dimension two, 



£ (2) = - j V" - ^A") 2 + *f(D„UD»rf + X f/t + x t|/) + C(Q R UQ L rf) , (1) 

where ( ) denotes the trace in flavor space, F^ u = d^A v — dyA^ the photon field strength 
tensor, A the gauge-fixing parameter (from here on, we work in the Lorentz gauge A = 1), 
-D M the generalized covariant derivative, 

D„U = d„U - ifa + a^ + QAJU + iU (v„ - a M + QA^) , (2) 

and Q the quark charge matrix, 

e ' 2 -°iH( 3t3 + 1 )' P) 

In what follows, we work in the so-called a-model gauge, 




U{x) = a{x) 1 + iT- n(x)/F , a(x) = y/l - tt 2 (x)/F 2 . (4) 

The third term in Eq.(|l]) is the standard non-linear a-model coupled to external sources 
(we neglect here singlet axial components and set (o M ) = 0) ||. In particular, we have 
X = 2i?o(s + ip) and the external scalar source contains the quark mass matrix, s(x) = 
diag(m u ,m^) + . . .. The last term in the dimension two Lagrangian is the lowest order 



chiral invariant term one can construct from pion and photon fields [0. To make it 
invariant under chiral SU(2)^xSU(2)^ transformations, we have introduced the spurions 
Ql,Ri which have the following transformation properties, 

Qi^giQigh gieSU(2) I} i = l,r. (5) 

At a later stage, one sets Ql = Qr = Q- The constant C can be calculated from the 
neutral to charged pion mass difference since this term leads to (6M 2 ) em = 2e 2 C/F 2 . This 
identification is based upon the fact that the quark mass difference rrid — m u only gives a 
tiny contribution to M 2 + — M 2 jlB] . It can therefore be expected that these quark mass 
effects are also tiny for the elastic tttt scattering amplitude. Consequently, we will work 
in the isospin limit = m u = rh in what follows. For later convenience, we introduce 
the dimensionless constant 

Z = C/F* = 0.89 , (6) 

with F = 88MeV ||]. As already stated before, we count the external vector and axial- 
vector fields as well as the charge matrices Q,Ql,Qr as 0(q) and the photon field as 
(9(1). This has the advantage of a consistent power counting between the strong and 
electromagnetic interactions, i.e. e ~ q and one has terms of dimension two, four and 
so on. Here, dimension two means either order q 2 or e 2 and similar at higher orders. 
The construction of the generating functional proceeds along standard lines making use 
of heat kernel techniques for elliptic Euclidean differential operators. We follow here the 
approach outlined in || and combine the expansion around the classical solutions for the 
meson and photons fields, U = U cl + iu^ u/F + . . . and A^ = A^ + e^, respectively, in one 
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i 


Oi 




1 


(QUQU^) 2 


3/2 + 3Z + 12Z 2 


2 


(QUQW) (D^UD^W) 


2Z 


3 


(UW^UQ) 2 + (D^UWQ) 2 


-3/4 


4 


(UW^UQ) (D»UWQ) 


-2Z 


5 


(IQ, [D^QjjWD^U) - ([Q, [D^QjjD^UW) 


-1/4 


6 


([D^Q}U[D^,Q}W) 





7 


(QUQW) (xf/t + Ux \) 


1/4 + 2Z 


8 


{{U^x - X ] U) (WQUQ - QWQU)) 


1/8 -Z 


9 


(Q 2 ) (QUQW) 


-3 - 3Z/5 - 12Z 2 /5 


10 


(Q 2 ) (D^UD^W) 


-27/20 - Z/5 


11 


(Q 2 ) (xt/ f + u x ] ) 


-1/4 - Z/5 


12 







13 


(Q 2 ) 2 


3/2- 12Z/5 + 84Z 2 /25 



Table 1: Counterterms and their /3-f unctions. 



set of fluctuations variables rj = . . . , £ 3 , e°, . . . , e 3 ). One then expands the generating 
functional 



exp[iZ(s,p,v l ^,a ljl )] = J\f J [cL4 M ] [dU] expji J d 4 x (C 2 + A) 



(7) 



up to second order in the fluctuations 77. Using dimensional regularization, the divergences 
can be extracted in a straightforward manner. For the two flavor case at hand, the em 
part of the dimension four counterterm Lagrangian takes the form 



13 



^em — *^kiOi , 



(8) 



i=i 



with the Oi monomials in the fields of dimension four. The low-energy constants fcj absorb 
the divergences in the standard manner, 



ki 
L 



KiL + kl(fj) 

,d-A r I 



(9) 
(10) 



with fj, the scale of dimensional regularization and d the number of space-time dimensions. 
The explicit expressions for the operators Oi and their /3-functions Ki are collected in 
table 1. The are the renormalized, finite and scale-dependent low-energy constants. 
These can be fixed by data or have to be estimated with the help of some model. 
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As it is the case for the hadronic LECs in the two-flavor case Q], one can introduce 
scale-independent couplings, called k^, via 



32tt 2 



L + In ■ 



Notice that we choose the neutral pion mass as the reference scale. This is natural since 
the neutral pion mass is almost entirely a hadronic effect, in contrast to the charged one. 
Some remarks concerning these results are in order. To arrive at the structures given in 
the table, we have made use of the equations of motions for the classical U field as well 
as of the Cayley-Hamilton matrix relations. This minimizes the number of independent 
terms. The last term in the table ~ (Q 2 ) 2 does not contain any pion field and is only 
needed for the complete renormalization, i.e. it does not influence physical processes. The 
same applies to Ou. Similarly, the operators 09,10,11 only lead to renormalizations of C, 
Fq and Mo, with Mq = 2rhBo the leading term in the quark mass expansion of the pion 
mass. A typical difference between SU(2) and SU(3) is the operator On, which has a 
vanishing /3-function in SU(3) ||. Note furthermore that symmetry does allow for terms 
which are odd in powers of Q, like e.g. a term ~ (Q 2 UQU^). Such terms are, however, 
unphysical since the electric charge always appears in even powers for physical processes. 
Consequently, we did not further study such terms. We also have repeated the SU(3) 
calculation and found that the /^-functions for the terms if 15,16,17 are incorrectly given in 
the existing literature. The correct ones read (in the notation of Urech 

E15 = 3/2 + 3Z + UZ 2 , £ 16 = -3 - 3Z/2 - Z 2 , S 17 = 3/2 - 3Z/2 + 5Z 2 . (12) 

This difference in the terms ~ Z 2 can be traced back to the fact that the operator 
a ab , which appears in the four-dimensional Euclidean one-loop functional, has to be 
symmetrized. This was not done in the part ~ C in ||, i.e. when going from Eq.(15) to 
Eq.(21) in that paper. 

3. Before evaluating the pertinent em corrections, we collect here some basic definitions 
concerning the elastic tttc scattering amplitude (in the threshold region). Consider the 
process 7i a (p a ) + n b {p h ) — > n c (p c ) + TT d (p d ), for pions with isospin 'a, b,c,d' and momenta 
p a ,b,c,d- The corresponding Mandelstam variables are s = (p a + Pb) 2 , t = (p a — p c ) 2 , 
u = (p a — Pd) 2 with s + t + u = M 2 a + M 2 b + M 2 C + M 2 d . The scattering amplitude can 
be expressed in terms of a single function, denoted A(s, t, u) ,F^| 

T cd;ab = u j 5 ab 5 cd + ^ ^ ^ gacfibd + ^ ^ ^ fiadfibc _ ^ 

The chiral expansion of A(s,t,u) takes the form 

A(s,t,u) =A^{s,t,u) + A^{s,t,u) + 0{q 6 ) , (14) 

where A^ is of order q m and the symbol 0(q 6 ) denotes terms like s 3 , s 2 t, s 2 e 2 , se 4 , .... 
The form of A^ 2 \s, t,u) was first given by Weinberg [[| and the next-to-leading or- 
der terms by Gasser and Leutwyler [|]. For convenience, one splits this contribution as 

# 6 Res Urech has kindly informed us that he has checked his coefficients £15,16,17 finding agreement 
with our results. 

# 7 This is correct in the isospin-conserving case. The more general form of the amplitude in case of 
isospin violation will be given in [p2[. 
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A^\s, t, u) = B(s, t, u) + C(s, t, u), where B(s, t, u) collects the unitarity corrections and 
C(s,t,u) the (real) tree and tadpole contributions. To evaluate the one-loop graphs, we 
introduce the modified pion propagator 



Af(£) 



ab 



v 2 - M 2 -5M 2 {l-5 3a )} 



5M 2 



2e 2 C 
F 2 o 



(15) 



with I the pion four-momentum and 'a, b' isospin indices. The form of Eq.([TJ) is due to 
the fact that in this gauge, the operator ~ C(QUQU>) only contributes to terms which 
are quadratic in pion fields. To evaluate the unitarity corrections, i.e. the diagrams which 
give rise to the imaginary part of the scattering amplitude, we introduce a generalization 
of the commonly used "bubble function" J, 



J(q 2 ,M 2 ,A) = 



16vr 2 



-In 



(7_ 



In 



<7_ 



1 



(16) 



with 



a 



<7_ 



AM 2 A [A 
1 + -(-+2 

s s \ s 



(17) 



/ 4M 2 /A 
s \ s 



1 



-2 



4M 2 



(Ja 



A 
M 2 



(7 



1 



-2 



with A a quantity of dimension [mass ] like e.g. 5M 2 . For A = 0, one recovers the 
standard form of J since then a = o~_ = a + = (1 — 4M 2 /s) 1 / 2 . All loop integral can 
be expressed in terms of J(q 2 , M 2 , A) and some polynoms (modulo logarithms). In terms 
of physical processes, we have five reaction channels 



(a) 
(d) 



7r°7r°- 

7T°7r + 



7r°7r° 



7r°7r- 



(b) 
(e) 



7T 7T 
7T+7T + 



7r°7r° 



71 TV 



7T 7T 



7T 7T 



(19) 



For comparison with the data, one decomposes T cd ' ,ab into amplitudes of definite total 
isospin (I = 0, 1, 2) and projects out partial-wave amplitudes 2/(s), 



1 - 4M 2 /s r 



2i 



exp{2H8{(8)+ir,t(8)]}-l 



(20) 



with s = 4:(M 2 + q 2 ) and q the pion momentum in the cm. system.^ 8 ] Furthermore, 
I denotes the total angular momentum of the two-pion system. The phase shifts Sf(s) 
are real and the inelasticities r/((s) set in at four-pion threshold. Below KK threshold, 
s ~ 1 GeV 2 , they are negligible and will be ignored in what follows. Near threshold, the 
partial-wave amplitudes take the form 



ReT l I (s)=q 2l {a I l+ qH I l +0(q±)} 



(21) 



The coefficients aj are called scattering lengths, the bj are the range parameters. In the 
following, we will concentrate on these quantities. 

# 8 This holds only for the equal mass case. The generalization to the unequal mass case is obvious. 
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4. We are now in the position to evaluate the em corrections to the elastic tttt scattering 
amplitude. To one-loop order, one has of course the standard strong interaction graphs, 
i.e. tree graphs at orders q 2 and g 4 as well as one-loop digrams at 0(q A ). In the a-model 
gauge, these can be calculated according to standard methods, the only difference being 
the modified pion propagator, Eq . fljjD . The corresponding Feynman diagrams are shown 
in fig. la,b,c,d. There are two tree graphs at order e 2 . Diagram a) contributes to the 
charged pion mass shift and b) vanishes for the reasons discussed above. That is also the 
reason why there is no one-loop graph like c) with the insertion on the four-pion vertex. 
The next seven graphs in fig.l are of order 0(e 2 q 2 ). From the irreducible photon loop 
diagrams (fig.le,f,g) only the first one is non-vanishing since one has no ^Att and 2747? 
vertices in this particular gauge (this point was also stressed by Gasser |I7j and Ecker (T^j ) . 



Clearly, graph h) are only gives rise to wave function renormalization and i) vanishes in 
dimensional regualarization. Finally, the counterterms with exactly one insertion from 
C^l are depicted in figs.lj,k. From these two, only graph lj gives rise to a genuine em 
correction to the 7T7T scattering amplitude. For charged pions, there are in addition photon 
exchange graphs as shown in fig. 2. At tree level 0(e 2 ) one has the one-photon exchange 
diagram. To one-loop order e 4 and e 2 q 2 , there are six topologically different graphs. 
While diagram 2a is of order e 2 q 2 , the others are of 0(e 4 ). In particular, there are three 
types of two-photon exchange graphs with one and two intermediate pion propagators, 
figs.2b,c,d, in order. The last two diagrams in fig.2 are simply vertex and self-energy 
corrections. As one easily convinces oneself, the photon exchange graphs like le, 2a, . . ., 
2f lead to IR divergences. These are cancelled in the standard fashion by considering 
the radiation of very soft photons in the initial and final states.* 9 For comparison with 
experiment, one can also use the standard Gamov factors [2~I] to remove the Coulomb 



enhancement in the initial and the final state for charged pions. What we are really after 
are the electromagnetic effects once these "kinematical" em effects are removed. Space 



forbids here to discuss these matters in detail and we refer to ref. f° r a comprehensive 
treatment. In what follows, we consider only the scattering of neutral pions, i.e. the 
reaction (a) in Eq.([L9[). 

5. For 7r°7r° — > 7r°7r°, we have s + t + u = 4M 2 and the leading order Weinberg term 



reads 



A^\s,t,u)= 5 -^ (22) 



in terms of the physical values M n o = 134.97 MeV and F w = 92.5 MeV.^ The respective 



shifts from the lowest order values M and F are accounted for in the next-to-leading 
order contribution C(s,t,u). It is important to stress that the lowest order amplitude 
is therefore not affected directly by the em corrections, only indirectly through the pion 



mass shift. This result is, of course, well-known [ l7| . The unitarity corrections take a 



form similar to the purely strong interaction results of 
1 



B(s, t, u) 



6Fi 



6(s - M 2 y J(s, M n+ ,0) - 3(s 2 - AsM 2 + 3M 4 ) J(s, M n o,0) 



# 9 A lucid discussion of such infrared effects is given in chapter 13 of ref.|T^]. A detailed study of 
radiative four-meson amplitudes in CHPT can be found in f20|| . 

# 10 Note that the extraction of this value for F„ includes one-loop radiative corrections p3|. 
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+ [t(t -u)- 2M 2 t + 4M 2 u - 2M* ]J(t, M n o, 0) 
+ [u(u - i) - 2M> + AM 2 t - 2M T 4 „]J(«, M^o, 0) 



(23) 



The next -to-leading order tree and tadpole contributions read 
1 



C(s, t, u) 



96vr 2 Fi 



2(/i -fa- 2Ml f + (l 2 - ^){s 2 + (t- uf] 



3M^ h 



+12M 2 {s - M 2 ) k - 12M> + 15M* 



1 



167T 2 Fi 



In 



7T~ 



M 2 + (3s-4M 2 ) + (s-M 2 



2 \2 



M\ 

+e 2 F, 2 (3s - 4M 2 



-h + 2Zk 4 



(24) 



The last two terms deserve some discussion. First, the term ~ ln(M 2 +/M 2 ) is due to 
the fact that we normalize the scale-independent SU(2) low-energy constants li to the 
neutral pion mass and that in the loops one has neutral as well as charged pion pairs 
propagating. The last term in Eq.([24|) is the novel em contribution from C^. To arrive 
at these results, we have used the low-energy expansions of the pion mass and the pion 
decay constant ,P TT 



Ml 



Ml 



1 - 



Ml 



F 1 + 



32tt 2 FI 
he 2 
72tt 2 

Ml 
16vr 2 F 2 

5e 2 

+ 1447T 2 



^3 + 



Ml 



2 '3 



■In 



Ml 



+ 



5e 



2 r 



16vr 2 F 2 Ml 72n 



Z 27\ 

2Z k 2 - ( - H )k w 

- 1 5 20/ 



32tt 2 



2Z+i) fc 7 



3 fc 3 + AZ h 



Z 



(25) 



I - 



Ml 



In 



Ml 



32tt 2 F 2 M 2 



^ 27 \ 

2Z k 2 - ( - H Uio 

' 1 5 20/ 



+ 



64tt 2 



3 A; 3 + 4Z fc 4 



(26) 



where M 2 = M 2 + 5M 2 , see Eq.(f[J). Notice that in B(s,t,u) and C(s,t,u) we have 
set F = F n , M = M n o and M + = M w + since these differences are of order q 2 and thus 
beyond the accuracy of the calculation presented here. As a check, we recover the result 
for + of |24| for one common pion mass and setting e = 0. From the amplitude 
given in Eqs.(22,23j24]), it is straightforward to calculate the pertinent scattering lengths 



a (00; 00) and effective ranges &o(00; 00) numerically or analytically [p2] . These are related 
to the ones in the isospin basis via 



a (00;00) 



a 2 , 6 (00;00) 



b°o + ~b 2 



o • 



(27) 



One can show analytically that for the process 7r 7r 
contribute. 



7r°7r°, the operators ~ hi do not 



# n In this letter, we identify the neutral with the charged pion decay constant. A more thorough 
discussion of the em effects on the AA correlators is given in . 



S 



For the numerical evaluation, we use the central values of the strong low-energy con- 



-1.7, 



6.0 together with 



2.9, 



4.3 and the values 



stants Ei from [^3J, t x - 
for F n , M n o and M n + given above. The corresponding S-wave threshold parameters are 
given in table 2 in comparison to the strong one-loop results and the experimental data. 





a (00;00) 


6 (00;00) 


e = 
e ^ 


0.0360 
0.0340 


0.0302 
0.0412 


Exp. 


0.056 ±0.027 26||27] 


0.029 ± 0.044 H 



Table 2: S-wave threshold parameters with em corrections at one loop (e ^ 0) compared 
to the hadronic one-loop results (e = 0) in units of the inverse neutral pion mass. The 
data should only be considered indicative since they mostly stem from processes involving 
charged pions. The uncertainties are added in quadrature. 

For the scattering length a® the effect of the em corrections is of the order of 5%, still 
approximately a factor of two smaller than the strong two-loop correction ||. For the 
range parameter bo, however, we observe an 36% increase. This is due to unitarity cusp 



at s = 4M^+, which is expected to scale as J M% + — M^ /M w + ~ 26%. The dominant 
isospin-violating effect is thus entirely given through the charged to neutral pion mass 
difference, similar to the case of the electric dipole amplitude in neutral pion photopro- 
duction off nucleons 



6. To summarize, we have constructed the generating functional for two-flavor chiral per- 
turbation theory including the effects of virtual photons in the one-loop approximation. 
Counting the electric charge as a small momentum, there are in total 13 terms contribut- 
ing to the em Lagrangian at next-to-leading order (some of these are only needed for 
renormalization). As an application, we have considered the em corrections to the elastic 
7T7T scattering amplitude and given numerical estimates for 7r°7r° — > 7r°7r°. The charged 
to neutral pion mass difference produces an pronounced effect on the S-wave effective 
range. In a forthcoming publication, we will present results also for the channels involv- 



ing charged pions [22] and include the effects due to the quark mass difference ~ — m v 
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Figures 




a) b) c) d) 




e) f) g) 




h) i) j) k) 



Figure 1: Graphs contributing to the em corrections to elastic tttt scattering. Solid 
and wiggly lines denote pions and photons, in order. The heavy dot and the box 
refer to insertions from the em counterterms at order e 2 and e 4 , respectively. Graphs 
b), f) and g) vanish in the a-model gauge, i) vanishes in dimensional regularization. 
Crossed graphs are not shown. 
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e) f) 



Figure 2: One and two-photon exchange graphs contributing to the em corrections 
for elastic scattering of charged pions at one loop. While diagram a) is of order e 2 q 2 , 
the others are of 0(e 4 ). Crossed graphs are not shown. For notations, see fig.l. 
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